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1.  Introduction 


Considerable  attention  1 1-19)  has  been  focused  on  the  photo-  or  radiation- 
induced  transient  hole  conductivity  in  amorphous  Si02  films  used  as  the 
gate  insulators  in  metal-oxide-semiconductor  (MOS)  structures.  A  prin¬ 
cipal  motivating  factor  for  this  work  was  the  realization  that  the  short-term 
postirradiation  transient  charge  and  current  response  of  MOS  systems  is 
dominated  by  hole  transport  through  the  Si02  films  16-10).  It  has  become 
apparent  that  this  hole  transport  is  rather  anomalous  in  nature.  In  par¬ 
ticular,  there  is  a  large  dispersion  in  the  transit  times  of  the  holes  through 
the  oxide.  Based  on  the  work  of  a  number  of  researchers,  the  general 
physical  picture  of  the  charge  response  of  Si02  MOS  systems  that  has 
emerged  is  the  following.  When  a  MOS  structure  is  exposed  to  ionizing 
radiation,  electron/hole  pairs  are  created  in  the  Si02  film.  Depending  upon 
the  applied  electric  field  across  the  Si02,  a  fraction  of  the  radiation¬ 
generated  carriers  undergoes  initial,  or  geminate,  recombination  |9, 
20-25  J.  This  fraction  is  dose  to  unity  at  low  fields  and  approaches  zero  at 
high  fields  (>3  MV/cm).  Those  electrons  escaping  geminate  recombina¬ 
tion  are  rapidly  swept  out  of  the  oxide,  in  times  on  the  order  of  pico¬ 
seconds  for  l(X)-nm  films  |26-29|.  Most  of  the  holes  remain  behind  near 
their  point  of  generation  and  cause  negative  flatband  voltage  or  threshold 
voltage  shift ,  in  the  MOS  structure.  Over  a  period  of  time  the  holes 
execute  a  relatively  slow  temperature-  and  field-activated  transport  |6-19, 
22 1  through  the  oxide  and  are  removed  at  an  interface  (at  the  Si  interface 
under  positive  bias  applied  to  the  gate  electrode).  This  hole  transport  is 
highly  dispersive  in  nature,  taking  place  over  many  decades  in  time,  and 
extending  in  the  tail  regions  to  seconds  at  room  temperature  and  much 
longer  at  lower  temperatures. 

The  origin  and  proper  description  of  the  dispersion  in  the  hole  transit  times 
has  been  a  topic  of  active  investigation.  In  our  initial  work  on  the  hole 
transport  properties  of  amorphous  Si02  |6.  1()|,  we  showed  that  the  time 
dispersion  in  the  time  regime  from  1()-4  to  10-'  s  could  be  well  described  in 
terms  of  a  stochastic  model  of  transport  in  disordered  materials  based  on  a 
continuous-time  random  walk  (CTRW)  |30-36|.  In  the  original  formula¬ 
tion  of  this  model  131-33,  35|.  the  transport  was  assumed  to  occur  via 
hopping  between  localized  gap  states  distributed  randomly  in  space.  The 
dispersion  was  attributed  to  a  large  distribution  of  intersite  hopping  times 
due  to  the  fluctuations  in  the  distance  between  neighboring  sites.  In  later 
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work  on  the  temperature  and  field  dependencies  of  the  transport  [11,  15|, 
we  suggested  that  the  microscopic  charge  transfer  mechanism  was  polaron 
hopping  between  the  localized  states.  We  also  showed  that  the  charac¬ 
teristic  transit  time  scaled  superlinearly  with  oxide  thickness  in  accordance 
with  the  prediction  of  the  CTRW  model  [17|.  R.  C.  Hughes  |9,  1 4 1  also 
interpreted  his  data  on  the  radiation-induced  conductivity  of  Si02  films  at 
late  times  (>10-6  s)  in  terms  of  the  CTRW  model.  At  earlier  times  he 
concluded  that  the  transport  proceeded  via  intrinsic  band  polaron  hopping. 
On  the  other  hand,  Curtis  et  al  1 12,  13!  analyzed  similar  experiments  on 
radiation-induced  currents  in  Si02  in  terms  of  a  multiple-trapping  mode!  in 
which  the  carriers  move  in  band  states  between  trapping  events,  and  the 
dispersion  arises  from  a  broad  exponential  decaying  distribution  of  trap- 
state  energy  levels  relative  to  the  valence  band  edge. 

Later,  several  workers  [37,  381  showed  the  mathematical  equivalence  of 
CTRW  and  multiple  trapping  (MT)  models  of  transport.  Transport  via 
hopping,  trap-controlled  hopping,  or  multiple  trapping  can  all  be  included 
as  special  cases  within  a  generalized  CTRW  framework  which  assumes 
only  a  broad  distribution  of  event  times.  The  central  question,  then,  is  not 
MT  versus  CTRW  models  of  transport,  but  rather  concerns  the  micro¬ 
scopic  origin  of  the  dispersion.  That  is,  does  the  dispersion  arise  from 
fluctuating  energy  levels  (such  as  a  distribution  of  activation  energies  in 
conventional  multiple  trapping),  or  does  it  arise  from  fluctuations  in  the 
overlap  integrals  for  hopping  (as  would  occur  for  a  random  spatial  dis¬ 
tribution  of  hopping  sites)?  A  key  issue  pertaining  to  the  origin  of  the 
dispersion  is  that  of  the  universality  of  the  transport,  particularly  with 
respect  to  temperature  |36|.  In  essence,  if  the  amount  of  dispersion  is 
invariant  with  respect  to  variations  in  temperature,  then  the  transport  is 
said  to  exhibit  universality  with  respect  to  temperature.  In  this  case  the 
dispersion  arises  from  fluctuations  in  overlap  integrals  for  hopping  and  not 
from  a  distribution  of  activation  energies.  On  the  other  hand,  if  the 
amount  of  dispersion  changes  with  temperature,  the  transport  is  not 
universal  with  temperature  and  the  dispersion  arises  at  least  in  part  from  a 
distribution  of  activation  energies.  The  earlier  studies  were  not  really 
definitive  on  the  issue  of  universality,  due  in  part  to  a  fairly  significant 
amount  of  long-term  deep  trapping  of  holes  which  tends  to  mask  the 
time-dependent  transport  behavior.  Thus  it  is  of  interest  to  address  the 
issue  of  universality  in  "clean"  amorphous  SiO,  in  which  there  is  a  low 
concentration  of  deep  trapping  sites. 


Accordingly,  in  this  report  we  present  the  results  of  a  detailed  study  on  the 
temperature,  field,  and  thickness  dependencies  of  the  transport  in  one 
particularly  clean  set  of  MOS  oxides.  In  these  samples,  less  than  2  percent 
of  the  holes  are  trapped  permanently  (i.e.,  for  long  times  compared  with 
the  characteristic  transport  time)  within  the  oxide  at  deep  trapping  sites. 
Hence,  these  samples  are  ideal  for  basic  transport  studies.  The  measure¬ 
ments  show  that  in  all  cases-for  the  entire  range  of  temperature,  field,  and 
thickness  investigated— the  time  dispersion  of  the  transport  can  be  ade¬ 
quately  described  by  a  single  value  of  the  CTRW  disorder  parameter,  a  = 
0.25  ±  0.03.  This  is  the  central  result  of  the  present  study,  and  indicates 
rather  strict  universality  of  the  hole  transport  with  respect  to  temperature 
(T),  field  (£),  and  thickness  (L)  in  clean  Si02.  That  is,  changing  these 
parameters  simply  changes  the  time  scale  for  the  response  but  not  the 
shape  or  overall  dispersion  of  the  response  curves.  In  more  quantitative 
terms,  the  response  may  be  characterized  by  the  functional  form  f(t/x\ a) 
where  the  disorder  parameter  a  describes  the  shape  (or  dispersion)  of  / 
and  is  independent  of  temperature,  field,  and  thickness.  These  parameters 
enter  only  in  the  characteristic  time  scale  for  the  response,  x  =  x(T,E,L). 
The  universal  behavior,  particularly  with  respect  to  temperature,  has 
implications  for  the  nature  of  transport,  namely,  that  the  transport  proceeds 
via  hopping  between  localized  sites.  The  origin  of  the  dispersion  lies 
primarily  with  a  distribution  of  intersite  transfer  integrals  as  would  arise, 
for  example,  from  a  random  spatial  distribution  of  hopping  sites  or  from  a 
distribution  of  bond  angles  in  the  network.  We  also  summarize  the 
evidence  that  indicates  that  the  microscopic  process  is  that  of  polaron 
hopping  between  localized  sites. 

In  the  next  section  we  briefly  review  the  general  concepts  underlying 
dispersive  transport  and  discuss  the  applications  of  the  CTRW  formalism 
as  a  general  theoretical  framework  for  its  description.  In  section  3  we 
discuss  the  experimental  techniques  used  in  the  measurements  as  well  as 
provide  a  description  of  the  samples.  The  experimental  results  are 
presented  in  section  4,  categorized  according  to  detailed  characterizations 
of  the  temperature,  field,  and  thickness  dependencies  of  the  transport.  In 
section  5  the  results  are  analyzed  in  depth  in  terms  of  the  CTRW  model, 
with  special  emphasis  on  exhibiting  the  universal  character  of  the 
transport.  In  the  last  section  wc  summarize  and  discuss  the  general 
implications  of  the  results. 


2.  Dispersive  Transport:  Theoretical  Review 

In  this  section  we  briefly  review  the  underlying  concepts  of  dispersive 
(non-Gaussian)  transport  in  disordered  solids  that  are  pertinent  to  the 
present  study.  For  more  information  concerning  details  of  the  mathemati¬ 
cal  formalism  or  for  a  review  of  the  experimental  situation  in  a  number  of 
solids,  the  reader  is  referred  to  the  excellent  review  article  |36|  by  Pfister 
and  Scher  (PS).  We  will  discuss  the  application  of  the  CTRW  formalism 
to  dispersive  transport,  including  the  development  of  approximate  solu¬ 
tions  with  which  to  compare  quantitatively  the  experimental  results  of  the 
charge  response;  and  we  discuss  the  basic  microscopic  transport 
mechanisms  that  can  be  described  within  the  CTRW  framework  and  their 
implication  for  interpreting  the  experimental  results. 

2.1  Non-Gaussian  Transport 

The  major  feature  which  distinguishes  transient  transport  in  disordered 
solids  is  the  broad  distribution  of  individual  (microscopic)  event  times 
which  extend  into  the  time  range  characteristic  of  the  experimental  ob¬ 
servations,  i.e.,  into  the  time  range  necessary  for  the  fastest  carriers  to 
transit  through  the  sample  |35-36|.  That  is,  the  distribution  of  event  times 
extends  to  or  further  than  the  transit  time  tf  .  This  leads  to  some  unique 
transport  characteristics.  In  particular,  the  transport  is  distinctly  non- 
Gaussian.  For  conventional  transport  obeying  Gaussian  statistics,  the 
dispersion  (o2  =  {(x  -  (x))2))  and  the  mean  displacement  ((*})  of  a  carrier 
packet  obey  the  relations  o  «  tm  and  (x)  oc  and  hence  o/(x)  ~  Cm.  In 
contrast,  dispersive  transport  in  disordered  solids  indicates  that  the  time 
dependence  of  the  dispersion,  or  spread,  of  the  carrier  packet  is  the  same 
as  the  mean  displacement,  i.e.,  o/(x)  =  constant.  In  essence,  while  some 
carriers  transit  the  sample  very  rapidly  via  a  succession  of  rapid  events 
(e.g.,  hops),  other  carriers  are  immobilized  near  their  generation  point  for 
times  of  the  order  of  or  greater  than  the  transit  time  of  the  fastest  carriers. 
To  quote  PS  [36),  "Under  extreme  non-Gaussian  conditions,  the  distribu¬ 
tion  of  the  carrier  packet  grows  asymmetrically  featuring  a  leading  edge 
penetrating  deep  into  the  bulk,  while  the  maximum  of  the  charge  density 
moves  only  slowly  out  of  the  generation  region.  For  such  asymmetric 
carrier  propagation,  the  spread  and  the  mean  position  have  the  same 
dependence,  hence  o/(x)  =  const."  Such  broad  distributions  of  event  times 
can  be  envisioned  in  hopping  transport,  where  small  fluctuations  in  either 
the  intersite  hopping  distance  or  activation  energy  (or  even  in  bond  angles) 
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can  produce  large  variations  in  the  hopping  transfer  integrals  and  hence 
lead  to  large  variations  in  hopping  times.  Similarly,  for  the  case  of 
transport  mediated  by  traps  (multiple  trapping),  relatively  small  variations 
in  the  trap  energy  levels  can  lead  to  a  broad  distribution  of  release  times. 

The  constancy  of  o/(x)  leads  to  several  unique  experimental  manifest¬ 
ations  of  dispersive  transport.  First  is  the  universal  character  of  the 
transport  in  which  the  shape  of  charge-response  curves  (e.g.,  current 
traces)  is  independent  of  the  transit  time  when  plotted  in  units  of  the  transit 
time.  Hence,  changes  in  any  experimental  parameter  which  do  not  affect 
the  dispersion  simply  change  the  time  scale  for  the  response  and  not  the 
shape  of  the  response  curves.  As  we  shall  see  for  hole  transport  in  amor¬ 
phous  Si02,  the  response  is  essentially  universal  with  respect  to  tempera¬ 
ture,  field,  and  thickness. 

A  second  unique  feature  of  dispersive  transport  is  that  of  an  apparent 
time-dependent  mobility,  or  alternatively,  a  thickness-dependent  mobility. 
As  time  progresses,  an  increasing  fraction  of  the  carriers  are  immobilized 
for  longer  times.  For  example,  in  hopping  transport  the  carriers  ex¬ 
perience  on  the  average  more  difficult  hops  with  time,  resulting  in  longer 
waiting  times  between  hops.  In  conventional  transport  language,  tiie 
effective  mobility  decreases  with  time,  or,  the  further  the  carrier  packet 
travels,  the  lower  the  effective  mobility  becomes.  As  an  immediate 
consequence,  the  transit  time  of  a  carrier  through  the  sample  varies  super- 
linearly  with  thickness.  This  is  in  sharp  contrast  with  conventional 
transport  based  on  Gaussian  statistics,  for  which  the  drift  mobility  is  a 
well-defined,  intrinsic  quantity,  and  the  transit  time  varies  linearly  with 
thickness. 

2.2  CTRW  Formalism 

Scher  and  Montroll  (SM)  |35|  applied  the  CTRW  formalism  to  the  cal¬ 
culation  of  transient  current  traces  for  dispersive  transport.  The  CTRW 
model  (as  developed  earlier  by  Montroll  and  Weiss  |3()|)  describes  a 
walker  hopping  randomly  on  a  periodic  lattice,  with  the  steps  occurring  at 
random  time  intervals  determined  by  a  hopping  time  distribution  function 
y(f).  The  existing  model  can  simulate  a  time-dependent  transport  process, 
in  the  presence  of  an  electric  field  and  absorbing  boundary,  in  which  the 
carriers  move  through  a  disordered  medium  with  a  broad  distribution  of 
event  (hopping)  times. 
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For  the  special  case  of  an  algebraic  distribution  function  of  the  form 
\|/(0  SM  provides  a  complete  mathematical  description  of  the 

transport  properties  in  terms  of  a  single  parameter,  the  disorder  parameter 
a,  which  assumes  a  value  between  zero  and  unity  and  which  is  determined 
by  the  detailed  microscopic  transport  process.  With  this  power  law  form 
for  \y(r),  both  SM  and  PS  have  been  able  to  successfully  account  for  the 
anomalous  transport  properties  in  several  inorganic  and  organic  amor¬ 
phous  materials.  In  particular,  this  form  for  the  distribution  function 
accounts  for  the  long  tails  in  the  transient  photocurrent  observed  in  thin 
films  of  these  materials.  For  the  experimental  situation  in  which  a  sheet  of 
charge  is  created  near  one  electrode  (e.g.,  by  ultraviolet  excitation)  and 
undergoes  stochastic  transport  to  the  opposite  electrode  under  the 
influence  of  an  applied  field,  the  predicted  current  has  the  form 

f r°-n)  ,  /</,  . 

/(/)  ,  (I) 


where  tr  is  the  transit  time  of  the  fastest  carriers  across  the  film.  In  the  SM 
treatment,  tr  is  given  by 


h 


exp(A JkT) 


(2) 


where  L  is  the  sample  thickness,  A  is  an  activation  energy  for  the  transport, 
and  1(E)  is  the  average  displacement  per  hop  and  depends  on  the  applied 
bias.  For  /  <  tr  the  current  decay  is  attributed  to  the  slowing  down  of  the 
carrier  packet  in  the  bulk,  which  in  turn  is  due  to  the  decrease  in  the 
effective  average  mobility  with  time.  For  t  >  tr  the  rate  of  decay  of  the 
current  increases  because  of  carriers  being  removed  at  the  absorbing 
boundary  (electrodes).  Note  that  the  sum  of  the  slopes  in  the  early-  and 
late-time  regimes  is  -2;  this  constitutes  a  critical  test  of  the  applicability  of 
a  power  law  form  for  the  distribution  function.  Also  note  from  equation 
(2)  the  superlinear  power  law  dependence  of  the  transit  time  with  thick¬ 
ness,  t,  ~  lJ'n ,  attributed  as  discussed  above  to  the  average  velocity  of  the 
carriers  decreasing  with  time.  SM  obtained  the  complete  solution  to  the 
transport  equations  of  the  biased  CTRW  only  for  the  special  (but  analyti¬ 
cally  tractable)  case  of  a  =  1/2.  For  other  values  of  a,  the  exact  solutions 
are  known  only  in  the  asymptotic  regions,  or  can  be  obtained  only  by 
detailed  computer  calculations. 
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We  point  oul  that  the  field  dependence  explicitly  indicated  in  equation  (2) 
is  applicable  only  at  low  fields.  For  the  high  fields  typically  encountered 
in  the  gate  oxides  of  MOS  systems  ( F.  >  1  MV/cm)  it  appears  that  HE)  is 
saturated  at  a  constant  \alue.  and  the  Field  dependence  results  from  a  field 
lowering  of  the  potential  harrier  heights  for  hopping  in  the  forward  direc¬ 
tion  Consequently ,  the  actu  ation  energy  A  in  equation  (2)  is  field  de¬ 
pendent  i  see  sect  5i 

Approximate  Solutions  to  (  TRW  Model 

In  previous  work  |b.  Id.  II.  15  IK)  we  have  successfully  applied  the 
(TRW  model,  using  a  power  law  distribution  function,  to  the  interpreta¬ 
tion  of  the  transient  charge  relaxation  in  thin  MOS  gate  oxide  films  (20  to 
X00  nmi  of  amorphous  SiO.  following  exposure  to  pulsed  ionizing  radia¬ 
tion  In  the  course  of  that  work  we  developed  simple  approximate  and 
analytically  tractable  solutions  to  the  biased  OTRW  model  for  arbitrary  a 
which  contain  all  the  essential  features  of  the  exact  solution  |39|.  The 
motivation  tor  developing  approximate  solutions  for  arbitrary  a  was 
severalfold  First,  the  range  of  (/-values  encountered  was  0.15  to  0.3.  so 
that  the  exact  solutions  for  a  =  1/2  were  inapplicable.  Second,  the 
penetrating  ionizing  radiation  generated  carriers  uniformly  throughout  the 
SiO,  films,  and  as  a  consequence  the  transition  region  between  early-  and 
late  time  regimes  was  greatly  smeared  out.  Hence,  to  quantitatively 
compare  the  experimental  curves  with  model  predictions  required  analyti¬ 
cal  solutions  for  the  charge  response  valid  for  all  time.  Finally,  the 
principal  physical  observable  used  to  monitor  the  charge  response  of  the 
SiO,  MOS  systems  was  the  flatband  voltage.  Al  r/?( / 5.  which  is  a  measure 
of  the  first  moment  of  the  time  dependent  charge  density  relative  to  the 
metal  (gatel  electrode,  and  the  description  of  Al  r/?(f)  is  somewhat  more 
complicated  than  for  the  current  teq  tin. 

In  this  report  we  compare  the  experimental  results  on  the  charge  response 
with  model  calculations  based  on  the  approximate  solutions.  For  details 
on  the  development  of  the  approximate  solutions  and  their  application  to 
numerical  calculations,  the  reader  is  referred  to  earlier  work  |39|.  In 
essence  we  require  appropriately  chosen  trial  functions  to  satisfy  the 
scaling  properties  and  the  general  requirements  of  the  theory.  These 
general  requirements  are  manifested  by  conditions  imposed  on  the  first 
three  moments  of  the  Green's  function  describing  the  charge  density 
response  of  the  svstem.  In  addition,  the  Green's  function  is  chosen  to 


exhibit  the  proper  (known)  asymptotic  behavior  at  early  and  late  times, 
which  is  dictated  in  part  by  the  presence  of  completely  absorbing  bound¬ 
aries.  In  the  calculation  of  physically  observable  quantities,  our  trial 
function  yields  exact  results  for  a  =  1/2,  where  the  exact  solution  is 
known. 

Charge  Response  of  MOS  Systems 

In  the  analysis  of  our  charge  relaxation  data,  we  are  confronted  with  the 
case  in  which  an  incident  pulse  of  ionizing  radiation  (12-MeV  electrons) 
produces  a  uniform  distribution  of  electron-hole  pairs  in  the  Si02  insulator 
film.  The  electrons  are  rapidly  swept  to  and  collected  by  the  electrodes, 
leaving  the  holes  to  execute  a  much  slower  stochastic  transport.  The 
principal  quantity  of  experimental  interest  is  the  radiation-induced  flatband 
voltage  shift  A VrBU)  in  the  capacitance-voltage  (C-V)  characteristics 
before  and  after  the  radiation  pulse.  The  quantity  AV^  is  proportional  to 
the  time-dependent  first  moment  of  the  charge  distribution  relative  to  the 
metal  electrode.  The  radiation  pulse  is  taken  to  occur  at  t  =  0.  Let  x  =  0  be 
the  metal/insulator  interface  and  .v  =  L  be  the  insulator/semiconductor 
interface.  Let  P(x.t;x  ')  be  the  Green's  function  for  the  charge  distribution 
at  x  and  r,  given  that  a  single  charge  was  created  at  x  =  x '  at  time  t  =  0.  If 
poU ')  is  the  initial  charge  distribution,  the  charge  density  at  x  at  time  t  is 
given  by 

p(.v.r)  dx'  pn(x')P(x,t:x')  .  (3) 


In  terms  of  p(jc,r),  Al 'FBU)  is 


A  V/nU)  = 


~r~  fti dx  r  p(r-r) 

OX 


where  Cox  -  £0X/L  is  the  oxide  capacitance  and  eox  is  the  dielectric  constant 
of  the  insulator.  Note  that  A Vr[iU)  is  negative  for  positive  charge  induced 
in  the  insulator.  For  the  case  of  uniform  ionization,  p(>U)  =  QJL,  where 
Qo  is  the  total  charge  per  unit  area  created  in  the  insulator  film  at  t  =  0.  In 
this  case  the  initial  flatband  voltage  shift  immediately  following  the 
radiation  pulse  is 


AV  >/,(()')  =- 


=  -1.9  x  1  (Y*L2Df(E) 


(5) 


where  D  is  the  radiation  dose  in  rad(Si02)  delivered  to  the  sample,  L  is 
expressed  in  nanometers,  and  f(E)  is  the  field-dependent  fraction  of 
carriers  escaping  geminate  (initial)  recombination.  For  the  radiation 
source  and  the  field  range  employed  in  this  investigation  (1  to  6  MV/cm), 
/(£)  ranges  from  -0.8  at  1  MV/cm  to  essentially  unity  at  6  MV/cm 
122-25|.  The  linearity  between  Qo  and  dose  contained  in  equation  (5)  is 
valid  at  relatively  low  doses  (<200  krad)  for  which  space-charge  perturba¬ 
tions  can  be  neglected.  At  higher  dose  levels,  the  relationship  between 
A\'rB(()+),  dose,  and  field  becomes  nonlinear  and  much  more  complex.  In 
the  present  study  the  dose  levels  are  sufficiently  low  that  equation  (5)  is 
applicable. 


To  compare  the  experimental  results  with  model  calculations,  we  employ 
the  approximation  scheme  outlined  elsewhere  |39|  for  the  Green's  func¬ 
tion  P(x.t;x ') .  Here  we  note  only  that,  corresponding  to  equation  (1 )  for 
the  asymptotic  behavior  of  the  current,  there  is  a  "critical  test”  of  the 
model  as  applied  to  the  asvmptotic  behavior  of  the  relaxation  of  AV  rM)  in 
the  case  of  uniform  ionization.  In  particular,  at  early  time  (r— >0)  the 
model  predicts 


-“'»«> -&(j -*-'1*) 


(6a) 


and  at  late  time  (r— *») 


-AV',,/,  (r) 


Qo  fh 

C, 


ox  t 


(6b) 


where  6,  and  b2  are  constants.  If  AV  rc(0+)  is  the  initial  shift  in  the  C-V 
characteristics  immediately  after  the  radiation  pulse,  the  charge  relaxation 
at  early  times  is  predicted  to  have  the  form  |Al (0+)  -  AV',/,  (0  I  ~  r2",  as 
/->(),  whereas  at  long  times  A l  (!)  decays  to  zero  as  Al  ,„  (/)  ~  in. 
Hence,  the  corresponding  test  of  the  model  would  be  the  observation  that 
in  terms  of  log-log  plots,  the  flatband  voltage  relaxation  at  early  times 
should  yield  a  slope  of  2a,  whereas  at  late  times  the  slope  should  be  -a. 
We  have  found  this  test  to  be  well  satisfied  by  our  earlier  studies  1 10.  18| 
of  hole  transport  in  SiO,.  We  will  apply  this  test  to  the  recovery  data  of 
the  present  study  and  show  that  model  calculations  using  the  a-value 
obtained  from  the  asymptotic  analysis  agree  well  with  the  time-dependent 
response  data  for  well  over  10  decades  of  time. 
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2.5  Microscopic  Transport  Mechanisms 

In  their  review  paper  Pfister  and  Scher  |36|  discuss  several  microscopic 
j  transport  mechanisms  which  can  lead  to  dispersive  transport  and  which 

|  can  be  treated  within  the  framework  of  the  CTRW  stochastic  model.  They 

|  are  indicated  schematically  in  figure  1.  The  first  is  multiple  trapping  in 

\  which  the  carriers  move  in  extended  band  states  between  trapping  events 

|  by  localized  states  (fig.  la).  The  localized  trap  states  are  distributed  in 

j  energy  (e)  which  gives  rise  to  large  fluctuations  in  the  thermal  release  rates 

from  the  traps,  and  hence  results  in  dispersive  transport.  In  this  case  the 
amount  of  dispersion  is  temperature  dependent,  caused  by  a  spread  in 
release  rates  which  vary  with  trap  depth  via  thermal  activation  exp (-e/fc'T) 
’  factors.  The  magnitude  of  the  relative  spread  in  release  rates  increases 

|  with  decreasing  temperature.  In  particular,  for  an  exponential  distribution 

of  traps  with  energy,  the  power  law  form  V|/(t)  ~  r!l+,,)  is  obtained  for  the 
CTRW  waiting  time  distribution  function,  where  the  disorder  parameter  a 
is  linearly  proportional  to  the  temperature. 


Figure  1 .  Schematic  of  possible  (a) 
charge  transfer  mechanisms, 
indicating  in  (a)  trap-mediated 
valence  band  hole  conduction,  and 
in  (b)  hopping  transport  via  direct 
hole  tunneling  between  localized 
*  trap  sites  within  Si02  bandgap. 


£  (b)  £ 

<£  C 
o«  ox 


The  second  transport  mechanism  is  hopping  transport  in  which  the  carriers 
move  directly  between  localized  trap  sites  via  phonon-assisted  tunneling 
processes  (fig.  lb).  Although  the  transport  in  this  case  is  in  general 
temperature  activated,  reflecting  the  potential  barriers  through  which  the 
carriers  tunnel  (e.g.,  polaron  hopping),  the  dispersion  arises  from  a  dis¬ 
tribution  of  overlap  (or  transfer)  integrals  for  hopping,  due  either  to  a 
random  distribution  of  intersite  hopping  distances  (as  treated  explicitly  by 
Scher  and  Lax  |31-32|)  or  to  a  distribution  of  bond  angles  in  the  amor¬ 
phous  structure,  which  can  modify  the  overlap  of  directed  orbitals.  In 
either  case,  in  contrast  to  multiple  trapping,  the  dispersion  should  be  only 
weakly  dependent  upon  temperature,  if  not  temperature  independent. 
Hence,  a  test  as  to  whether  pure  hopping  is  involved  would  be  the  univer¬ 
sality  of  the  response  curves  with  respect  to  temperature.  Of  course,  the 
dispersion  can  arise  from  a  combination  of  fluctuating  hopping-overlap 
integrals  and  activation  energies,  in  which  case  the  temperature  depend- 
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ence  of  the  dispersion  would  lie  somewhere  between  that  for  pure  hopping 
and  that  for  multiple  trapping.  In  this  case  the  stochastic  features  of  the 
transport  can  be  described  within  the  generalized  CTRW  framework,  but  it 
would  be  difficult  to  say  anything  definite  about  the  nature  of  the  micro¬ 
scopic  transport  mechanism. 

We  have  discussed  hopping  transport  in  a  rather  generic  fashion,  simply 
indicating  that  it  generally  involves  phonon -assisted  tunneling  transitions 
between  two  localized  trap  sites.  Now  we  focus  on  a  single  hopping 
mechanism-that  of  small  polaron  hopping-because  as  we  shall  see,  some 
of  the  specific  characteristics  (temperature  and  field  dependencies)  of  hole 
transport  in  Si02  are  best  explained  by  this  process.  A  detailed  discussion 
of  the  theory  of  polaron  formation  and  transport  is  beyond  the  scope  of  the 
present  work;  for  details  we  refer  the  reader  to  the  extensive  literature  on 
the  subject  (see  for  instance  references  40-43).  We  will,  however,  indicate 
the  essence  of  the  process  in  a  qualitative,  schematic  fashion,  and  briefly 
reiterate  some  of  the  main  experimental  manifestations  of  the  process. 

The  driving  mechanism  responsible  for  polaron  formation  is  a  strong 
interaction  between  a  charge  carrier  and  the  lattice,  such  that  the  total 
energy  of  the  system  is  lowered  by  a  distortion  of  the  lattice  (or  atom 
network  for  a  disordered  solid)  in  the  vicinity  of  the  carrier.  In  essence, 
the  carrier  polarizes  the  surrounding  medium,  and  this  polarization  then 
interacts  back  on  the  carrier.  If  the  interaction  is  sufficiently  strong  with  a 
large  distortion  of  the  lattice  in  the  immediate  vicinity  of  the  carrier,  then 
the  carrier  becomes  localized  at  a  particular  site.  This  is  referred  to  as  a 
small  polaron,  or  sometimes  also  as  self-trapping  of  the  carrier.  The 
localization  of  the  carrier  can,  in  general,  be  either  at  a  lattice  site  or  near  a 
defect  (trap  site).  Particularly  in  a  disordered  solid,  the  localization  may 
be  associated  with  bandtail  states  originating  from  the  disorder.  Such  a 
strong  carrier-lattice  interaction  leading  to  small  polaron  formation  is  a 
common  feature  in  polar  materials,  or,  more  generally,  in  materials  with  a 
large  polar  component  in  their  electronic  structures. 

In  the  simplest  small  polaron  theory  (see,  for  example,  Mott  and  Davis 
1 4 1  ]),  it  is  assumed  that  the  electronic  energy  varies  linearly  with  the 
configuration  coordinates  of  the  system  (relative  to  the  equilibrium 
configuration  in  the  absence  of  the  carrier),  and  the  elastic  strain  energy 
varies  quadratically.  Then,  the  net  reduction  of  the  total  system  energy 
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due  to  the  polarization  distortion  of  the  medium  is  -IV  ,  which  is  the 
resultant  combination  of  a  -2 Wp  reduction  in  the  electronic  energy  and  an 
increase  of  Wp  in  the  elastic  strain  energy. 

The  mechanism  of  charge  transport  via  polaron  hopping  is  depicted 
schematically  in  figure  2.  In  essence,  when  an  initially  empty  localized 
trap  site  (fig.  2a)  captures  a  carrier  (in  our  case,  a  hole),  the  total  energy  of 
the  system  is  lowered  by  a  distortion  of  the  lattice  around  the  trap  site  (fig. 
2b).  The  hole  "digs"  a  potential  well  for  itself,  i.e.,  is  self-trapped.  When 
the  hole  moves  through  the  insulator,  it  carries  with  it  the  potential  well 
arising  from  the  distortion  of  the  lattice.  The  transition  of  the  hole  be¬ 
tween  two  nearby  sites  occurs  via  the  intermediate,  thermally  activated 
state  shown  in  figure  2c,  for  which  thermal  fluctuations  of  the  system 
momentarily  bring  the  electronic  energy  levels  on  the  two  sites  into  coinci¬ 
dence,  and  the  hole  tunnels  from  site  (1)  to  site  (2).  The  final  state  of  the 
system  is  depicted  in  figure  2d,  with  the  hole  now  residing  in  site  (2). 


Figure  2.  Schematic  illustration  of 
polaron  hopping  between  two  nearby 
localized  trap  sites,  (1)  and  (2). 

(a)  Both  sites  unoccupied,  (b)  Hole 
localized  on  site  ( I ).  (c)  Direct 
tunneling  of  hole  from  ( 1 )  to  (2)  when 
energy  levels  arc  brought  into  coin¬ 
cidence  by  thermal  fluctuations. 

(d)  Final  state  with  hole  localized  on 
site  (2). 
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The  transition  probability  for  the  process  is  the  product  essentially  of  two 
factors.  One  is  the  spatial  overlap  integral  of  the  wave  functions  on  the 
two  sites  which  governs  the  tunneling  transition  between  the  two  sites.  It 
is  of  the  form  exp(-2|ki2),  where  xn  is  the  intersite  separation  distance 
and  (5  is  related  to  the  potential  barrier  between  the  sites.  The  second 
factor  is  the  probability  of  creating  the  intermediate,  activated  state  of 
coincident  energy  levels  shown  in  figure  2c. 

At  high  temperatures  (T  >  TDN ,  where  To  is  the  Debye  temperature  and  v 
is  a  numerical  factor  on  the  order  of  2  to  4),  the  energy  to  form  this 
coincident  state  is  supplied  by  the  thermal  fluctuations  (phonons)  of  the 
system.  Hence,  the  process  is  thermally  activated  with  an  activation 
energy  A  =  ( Wp  +  Wd)/2  where  Wp  is  the  net  reduction  in  the  system 
energy  associated  with  the  polarization  of  a  single  site,  and  Wd  is  a  disor¬ 
der  energy  associated  with  the  spread  of  the  intrinsic  energy  levels  of  the 
hopping  sites  (41,  44|.  We  note  that  the  effect  of  an  applied  field  is  to 
effectively  reduce  the  disorder  energy  by  an  amount  equal  to  the  average 
potential  drop  between  neighboring  hopping  sites. 

For  low  temperatures  ( T  <  T^/v),  the  activation  energy  associated  with  the 
transition  depicted  in  figure  2  decreases  from  (Wp  +  Wd)/ 2  to  the  (usually) 
much  smaller  value  Wd.  The  reason  for  this  is  that  the  coincidence  con¬ 
figuration  of  figure  2c  results  in  this  case  primarily  from  the  zero-point 
vibrations  of  the  system.  Basically,  the  total  energy  EL  associated  with  the 
lattice  vibrations  (phonons)  can  be  written  as  EL  =  E7Jt  +  Er,  where  E7J,  is 
the  quantum  mechanical  zero-point  energy  of  the  lattice  (the  lattice 
vibration  energy  at  zero  temperature)  and  Ef  is  the  energy  associated  with 
the  thermally  induced  fluctuations  of  the  system.  Since  T-tp  is  fixed, there 
is  a  crossover  point  in  temperature  at  which  E7J>  =  Er  that  occurs  at  a 
temperature  ~T0/v.  The  actual  temperature  region  of  the  transition  (i.e., 
the  value  of  v)  depends  primarily  on  details  of  the  phonon  spectra  of  the 
solid.  Above  the  transition  temperature  the  lattice  vibrations  are  primarily 
thermally  induced,  whereas  well  below  this  temperature  the  lattice  vibra¬ 
tions  are  primarily  due  to  the  zero-point  oscillations.  For  a  nonzero  disor¬ 
der  energy,  however,  some  thermal  component  in  the  activation  remains 
even  at  very  low  temperature,  because  some  net  transfer  of  energy  (i.e., 
absorption  or  emission  of  a  single  phonon)  between  the  lattice  and  the 
electronic  energy  is  required  for  the  transport  process  to  occur;  and  this  net 
energy  transfer  cannot  involve  E7j,  as  this  energy  component  is  a  constant. 


The  transition  from  a  thermally  activated  process  to  one  that  is  not  ther¬ 
mally  activated,  or  to  one  with  a  much  lower  activation  energy,  as  the 
temperature  is  lowered  below  some  value  (-T^/v)  is  a  classic  signature  of 
polaron  hopping  1 40-45 1.  As  we  see  in  section  4,  the  temperature  depend¬ 
ence  of  the  hole  transport  in  Si02  exhibits  such  a  transition.  Also,  the 
main  effect  of  an  applied  electric  field  is  to  reduce  the  activation  energy  of 
the  transport  in  a  manner  expected  in  a  polaron  hopping  process.  Further¬ 
more,  that  polaron  hopping  may  be  the  primary  microscopic  charge 
transfer  process  is  not  unexpected,  since  Si02  has  a  large  polar  component 
in  its  electronic  structure  associated  with  a  narrow  band  of  states  lying  at 
the  top  of  the  valence  band,  which  originate  from  nonbonding  2p  orbitals 
on  the  oxygen  atoms  |46-50|.  These  states,  which  are  not  part  of  the 
covalent  bonding  of  the  Si02  network  (the  orbitals  are  directed  away  from 
the  covalent  bonds),  are  easily  deformable  and  give  rise  to  a  large 
polarizability  of  the  structure. 
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IX 


3.  Samples  and  Experimental  Techniques 


Thin  film  SiO,  capacitors,  grown  as  a  part  of  a  radiation-hardening  effort, 
were  supplied  by  Hughes  Aircraft  Corporation.  The  samples  consisted  of 
wet-grown  SiO,  layers  deposited  on  (1(H))  n-  and  p-tvpe  Si  substrates  with 
vapor-deposited  A1  gate  electrodes.  The  wet-process  SiO,  was  grown 
with  pyrogenic  H,0  at  950°C  to  a  96.5-nm  thickness  and  annealed  in  N,  at 
925  C  for  20  minutes.  These  samples  are  especially  suitable  for  transport 
studies  because,  as  shown  by  late-time  measurements,  less  than  2  percent 
of  the  holes  are  permanently  trapped  in  the  oxide. 

For  the  thickness  dependence  studies,  two  sets  of  MOS  capacitors  w'ere 
used,  each  grown  on  n-type  Si  substrates  with  SiO,  gate  insulators  of 
various  thicknesses.  Both  sets  of  insulators  were  thermally  grown  in  a 
pyrogenic  H20  ambient  at  925*  0  and  annealed  in  dry  nitrogen.  For  the 
"etched-back"  set.  gate  oxidation  was  carried  out  simultaneously  on  all 
wafers  to  produce  nominal  1(H)- nm  films.  Some  of  the  films  were  then 
etched  to  yield  samples  of  nominal  70-.  60-.  and  40-nm  Si02  thicknesses. 
Since  the  initial  gate  oxidation  for  these  samples  was  the  same,  the 
etched-back  samples  were  expected  to  have  similar  bulk  oxide  and  Si02/Si 
interface  properties.  The  "as-grown"  sample  set  was  produced  by  varia¬ 
tion  of  the  gate  oxidation  time,  while  a  constant  oxidation  temperature  and 
ambient  were  maintained,  producing  samples  with  nominal  1(H)-,  80-,  60-. 
40- .  and  20-nm  SiO,  thicknesses.  These  samples  more  closely  represent  a 
practical  commercial  product  than  do  the  etched-back  samples. 

The  experiments  were  carried  out  by  using  the  electron  linear  accelerator 
(LIN AC)  at  the  Armed  Forces  Radiobiologv  Research  Institute.  The 
L1NAC  pnxluced  a  nominal  15-MeV  electron  beam  with  a  pulse  width  of 
4  ps.  The  dose  delivered  to  the  samples  was  kept  below  20  krad(Si02)  in 
all  cases  to  minimize  perturbation  of  the  electric  field  in  the  SiO,  by  the 
radiation-produced  space  charge.  The  samples  were  irradiated  in  an 
evacuated  sample  holder  which  had  provisions  for  electron  beam 
dosimetry  and  for  control  of  sample  temperature  from  80  to  5(H)  K. 
Charge  relaxation,  or  transport,  was  measured  by  monitoring  the  shift  in 
the  Hatband  voltage  in  the  capacitance-voltage  (C-V)  characteristic 
between  the  preirradiation  value  and  the  value  at  logarithmically  increas¬ 
ing  time  intervals  after  the  radiation  pulse  from  0.2  ms  to  8(H)  s.  using  a 
fast  high  frequents  C-V  measurement  apparatus.  As  discussed  in  the 
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previous  section  (see,  e.g.,  eq  (4)),  the  flatband  voltage  shift  A VrB(t)  is 
proportional  to  the  time-dependent  first  spatial  moment  (relative  to  the 
metal-oxide  interface)  of  the  radiation-induced  charge  distribution  in  the 
oxide.  As  the  charge  transports  and  is  collected  at  the  elect  xles,  the 
flatband  voltage  relaxes  toward  its  preirradiation  value. 

A  block  diagram  of  the  fast  C-V  apparatus  is  shown  in  figure  3  1 17).  The 
MOS  capacitor  in  the  sample  holder  is  normally  maintained  under  a  preset 
dc  bias  during  and  after  a  radiation  pulse.  This  bias  voltage  is  set  by  a 
voltage  reference  and  routed  to  the  sample  through  a  series  of  electronic 
switches,  a  power  amplifier,  and  a  coupling  network.  The  MOS  sample 
capacitance  is  continuously  monitored:  a  low  level  of  rf  excitation  (typi¬ 
cally  at  5  MHz)  is  applied  to  the  sample  and  a  phase-sensitive  detector 
measures  the  reactive  component  of  the  resultant  current  through  the 
sample  capacitance.  Measurements  of  AVFB  as  a  function  of  time  after  a 
radiation  pulse  are  controlled  by  the  timing  generator.  When  activated  by 
a  signal  from  the  radiation  source,  this  unit  produces  a  series  of  trigger 
pulses  at  logarithmically  increasing  time  intervals  (in  a  2-4-8-20  sequence) 
typically  from  0.2  ms  to  800  s  after  the  radiation  pulse.  Upon  receipt  of  a 
trigger  pulse,  a  \'FB  measurement  is  performed  by  the  \'FB  digitizer  subsys¬ 
tem  in  the  following  sequence:  (1)  The  dc  bias  on  the  sample  is  inter¬ 
rupted  and  the  accumulation  capacitance,  Cox,  of  the  sample  is  briefly 
sampled.  (2)  A  voltage  ramp  is  applied  to  the  sample  and  passage  of  the 
sample  capacitance  through  a  preselected  value  of  CFB/C()X  (where  CFB  is 
the  flatband  capacitance)  is  detected.  (3)  The  corresponding  ramp  voltage, 
\'FB,  is  captured,  digitized,  and  stored  in  a  fast  memory.  For  the  experi¬ 
ments  reported  here,  a  complete  VFB  measurement  cycle  required  -100  ps. 
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4.  Results 


The  experimental  results  may  be  grouped  for  purposes  of  presentation  into 
three  general  categories:  (a)  temperature  characterizations,  where  the  field 
is  held  fixed  and  the  temperature  varied,  nominally  between  79  K  and 
room  temperature;  (b)  field  characterizations,  where  the  temperature  is 
held  constant  and  the  field  varied,  and  (c)  thickness  dependence  data, 
where  the  results  are  presented  for  oxides  of  several  thicknesses  at  the 
same  field  and  temperature.  The  results  for  (a)  and  (b)  are  presented  in 
detail  for  one  oxide  thickness  only. 

4.1  Temperature  Characterizations 

The  fiatband  voltage  recovery  data  are  shown  in  figures  4  to  7  for  oxide 
fields  of  1,  3.  5,  and  -1  MV/cm,  respectively.  In  each  case  the  data  for  a 
series  of  temperatures  in  the  nominal  range  between  liquid  nitrogen  and 
room  temperature  are  plotted  versus  log  time  from  0.2  ms  to  800  s  follow¬ 
ing  the  radiation  pulse.  The  data  are  normalized  to  the  initial  shifts, 
AVrfl(0+),  immediately  after  the  pulse  in  order  to  simplify  the  analysis  and 
comparisons  of  the  data.  The  actual  values  of  &VrBi.0+)  vary  with  dose 
and  field,  as  described  in  section  2.  The  data  are  plotted  in  the  negative 
direction  because  the  shifts  are  negative,  indicative  of  net  positive  charge 
induced  in  the  oxide  layer. 
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Figure  4.  Platband  voltage  recovery  data  following  Figure  5.  Normalized  Hatband  voltage  recovery  data 

pulsed  LIN  AC'  electron  irradiation  ol  46.5-nm  lor  oxide  field  of  f  MV/cm  for  temperatures  between 

wet-grown  oxide  capacitor  under  I  -MV/cm  oxide  X2  and  2lfS  K. 

field  for  scries  of  temperatures  between  1 24  and 

2^1  K  The  data  are  normalized  to  the  initial  Hatband 

voltage  shilv  immediately  alter  the  radiation  pulse. 


.^jV'jTrjir'  iT.rsr. 


iQo'  ood 


(T  246K) 

_  o 


□  08! 

□  O 

o  A  o  O 
Q  a  $  o 
aA«0 


(219K) A 
A 
A 


□ 

o°  > 

O  ° 


□  u  A  V  o 
□  Q  A  Q  O 
o  A0o 
a  a  0  o 

A  A  A  O 


□  __  D 
Q  Dht8k, 

093*,°  Q  Q  , 


l  (161  K.)  q 

°° 

Q  Q  [140K]  Q  ° 


Q%  °° 


o0°  0 

V  (123KI/S 


O  v  { 1 23K)  Q 


Q°  o°°  XOO  00?  ° 


60”  v  O 

Cl  O  o  q  o  .  (109KI 


,0-»  10-3  10-2  10-1 


103  ,0  3  V  ,^-3V  V  ,0  ='  W7'  10° 


10  10  10 


TIME  AFTER  PULSE  <s| 


TIME  AFTER  PULSE  <s| 


Figure  6.  Normalized  Hatband  voltage  recovery  data 
for  oxide  field  of  5  MV/cm  for  temperatures  between 
79  and  250  K. 


Figure  7.  Normalized  Hatband  voltage  recovery  data 
for  negative  gate  oxide  Held  of -  I  MV/cm  for 
temperatures  between  109  and  246  K. 


There  are  several  points  to  note  about  the  data.  First,  it  is  very  evident 
from  all  the  data  sets  that  the  recovery,  or  hole  transport,  is  very  dispersive 
in  time,  taking  place  over  at  least  8  decades  in  time  (more  like  10  decades 
from  start  to  finish,  as  will  become  more  apparent  in  the  analysis  section). 
It  is  also  evident  that  the  recovery  is  strongly  temperature  activated,  with 
little  recovery  (at  least  for  oxide  fields  Eox  <  3  MV/cm)  occurring  at  the 
low  temperatures,  but  with  increasing  recovery  taking  place  as  the  tem¬ 
perature  increases,  and  with  the  recovery  essentially  complete  at  the 
highest  temperatures  (room  temperature  for  Enx  <  3  MV/cm,  250  K  at 
5  MV/cm).  The  almost  complete  recovery  at  late  times  and  high  tempera¬ 
ture  indicates  very  little  long-term  (permanent)  trapping  (<2  percent)  of 
the  holes  in  these  samples.  At  Eox  =  5  MV/cm  (fig.  7),  there  are  evidently 
two  distinct  temperature  regimes.  Above  about  140  K,  the  transport  is 
strongly  temperature  activated,  but  it  is  only  weakly  temperature  depend 
ent  below  140  K.  But  note  the  universal  appearance  of  the  data  in  all  cases 
(including  the  two  temperature  regimes  for  the  5-MV/cm  data):  namely, 
changing  the  temperature  affects  the  time  scale  of  the  transport  without 
significantly  changing  the  shape  of  the  recovery  curve.  These  points  are 
elaborated  in  greater  depth  in  section  5. 


The  data  in  figure  7  for  the  negative  bias  case  (-1  MV/cm)  are  similar  in 
appearance  to  the  +1  MV/cm  recovery  data.  There  are  differences, 
however,  specifically  with  respect  to  an  increased  speed  of  the  flatband 


voltage  recovery  (due  to  the  motion  of  holes  in  this  case  toward  the  metal 
gate  electrode,  relative  to  which  the  spatial  moment  is  taken  in  eq  (4)  for 
&VFB)  and  to  a  difference  in  the  shape  of  the  recovery  curves  at  early  time. 
These  points  are  more  easily  addressed  when  we  discuss  the  universal 
nature  of  the  transport  in  detail  in  section  5.  For  the  moment  we  note  that 
the  initial  flatband  voltage  shift  at  -1  MV/cm  is  also  negative  and  within 
~5  percent  of  the  initial  shift  for  1  MV/cm  normalized  to  the  same  dose. 
Since  electrons  are  very  mobile  in  Si02  and  are  swept  out  under  either  bias 
on  the  order  of  picoseconds,  this  result  is  consistent  with  uniform  genera¬ 
tion  of  the  holes  throughout  the  Si02.  At  low  temperature  and  field,  the 
holes  remain  frozen  very  near  their  generation  point  for  the  time  scale  of 
the  experiment  (800  s).  At  higher  fields  (5  MV/cm)  significant  recovery 
does  occur  during  the  experimental  observation  time  even  at  79  K,  which 
leads  now  to  the  field  dependence  of  the  transport. 

4.2  Field  Dependence  of  the  Recovery 

First,  we  present  data  taken  from  Figures  4  to  7  in  a  way  to  bring  out  the 
field  dependence  of  the  transport.  In  figure  8  we  show  three  sets  of 
recovery  data  with  each  set  showing  the  recovery  at  oxide  fields  of  1 ,  3, 
and  5  MV/cm  at  a  constant  temperature:  T  =  124  K  in  figure  8a,  165  K  in 
(b),  and  250  K  in  (c).  Again  the  speed-up  of  the  transport  with  tempera¬ 
ture  is  evident:  at  124  K  the  early  part  of  the  recovery  is  observed, 
whereas  at  250  K  the  latter  portion  of  the  recovery  is  observed,  and  the 
recovery  at  165  K  encompasses  the  mid  range  between  the  early  and  late 
stages.  The  field  dependence  of  the  recovery  is  evident  in  each  case.  But 
the  main  point  to  be  made  here  is  that  the  transport  becomes  more  field 
activated  as  the  temperature  is  lowered.  Note  that  the  total  spread  in  the 
recovery  between  1  and  5  MV/cm  increases  from  about  3  decades  at 
250  K  to  about  6  decades  at  124  K.  This  is  the  behavior  expected  for 
activated  transport  where  the  main  effect  of  the  field  is  to  reduce  the 
activation  energy  for  the  transport. 

In  figure  9  the  field-dependent  data  are  shown  for  T  =  79  K  for  positive 
fields  from  3  to  6  MV/cm,  for  which  case  the  transport  is  highly  field 
activated  1 15).  Below  3  MV/cm  there  is  essentially  no  transport  occurring 
in  the  measurement  time.  At  3  MV/cm,  some  transport  and  recovery  begin 
to  take  place  on  the  time  scale  of  the  observations,  and  as  the  field  is 
increased  above  3  MV/cm,  the  charge  transport  increases  strongly  until,  at 
6  MV/cm,  there  is  almost  complete  recovery  of  Al7rfl  at  800  s,  indicating 
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Figure  8.  Normalized 
Platband  voltage  recovery 
data  for  several  values  of 
oxide  Held  (1,1,  and  5 
MV/cm)  at  several 
temperatures:  (a)  at 
124  K;  (b)  at  165  K;  and 
(c)  at  250  K. 
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that  most  of  the  positive  charge  has  been  collected  at  the  electrode. 
Doubling  the  field  from  3  to  6  MV/cm  results  in  an  increase  in  transport 
speed  of  6  orders  of  magnitude.  Note  again  the  universal  appearance  of 
the  data:  the  main  effect  of  changing  the  field  is  to  shift  the  time  scale  for 
the  recovery  without  appreciably  affecting  the  shape--or  amount  of  overall 
dispersion  -of  the  recovery  curves. 


Figure  9.  Normalized 
flaiba-  j  voltage  recovery 
dau  .ollowing  pulsed 
LINAC  electron-beam 
exposure  for  96.5-nm 
oxide  at  80  K  and  for 
oxide  fields  from  3  to 
6  MV/cm. 
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Thickness  Dependence 

Figure  10  presents  high  field  (5-MV/cm)  transport  data  1 1 7 1  at  79  K  for 
three  oxide  thicknesses:  99,  72,  and  41  nm.  The  samples  in  this  case 
were  all  grown  to  the  same  thickness  (99  nm)  and  then  portions  of  the 
wafer  preferentially  etched  back  to  obtain  the  small  oxide  thicknesses. 
Data  were  also  obtained  for  a  59-nm  oxide  sample,  but  its  recovery 
tracked  so  closely  (though  shifted  slightly  to  earlier  time)  to  the  72-nm 
recovery  curve  that  it  was  omitted  from  the  figure  for  clarity. 


The  large  electric  field  for  the  measurements  in  Figure  10  causes  relatively 
rapid  annealing  of  &VrB  even  at  low  temperature  (80  K).  It  is  evident  that 
recovery  takes  place  more  rapidly  in  the  thinner  oxide  samples; 
however,  the  shapes  of  the  annealing  curves  do  not  differ  significantly. 
Therefore,  the  primary  effect  of  varying  oxide  thickness  is  again  a  transla¬ 
tion  of  the  recovery  curves  to  earlier  time  with  decreasing  thickness;  i.e., 
only  the  time  scale  of  the  hole  transport  is  affected  as  the  oxide  is  thinned. 
Thus,  this  "universal"  behavior  has  been  obtained  now  for  variations  in 
temperature,  field,  and  oxide  thickness.  We  illustrate  this  universality 
feature  more  graphically  in  the  analysis  section. 


The  data  in  figure  10  also  show  the  strong  superlinear  dependence  of  hole 
transit  time  with  thickness.  For  the  99-nm  sample,  AV^  anneals  to  one 
half  the  initial  flatband  shift  in  0.5  s.  For  the  41-nm  sample,  the  half¬ 
anneal  point  is  reached  in  0.015  s.  Thus,  a  reduction  in  oxide  thickness  of 
a  factor  of  2.4  shortens  the  recovery  time  by  a  factor  of  33.  Keep  in  mind 
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Figure  10.  Normalized 
Hatband  voltage  recovery 
data  following  pulsed 
LIN.  ■T'.  electron  beam 
irradiation  for  three 
thicknesses  (4 1 , 72,  and 
99  nm)  of  elched-back  dry 
oxide  samples  at  80  K  and 
applied  oxide  field  of 
5  MV/cm. 


that  for  conventional  charge  transport  models  based  on  a  constant 
mobility,  the  transit  time,  or  recovery  time,  would  simply  vary  linearly 
with  thickness. 

Figure  1 1  presents  normalized  &VpB  results  for  as-grown  Si02  samples 
with  thicknesses  of  96,  56,  38,  and  22  nm  irradiated  at  2(K)  K  with  an 
applied  field  across  the  oxides  of  1  MV/cm.  Again,  it  is  apparent  that  the 
AV FB  recovery  curves  shift  superlinearly  to  earlier  times  with  decreasing 
oxide  thickness.  In  this  case,  a  reduction  in  oxide  thickness  of  a  factor  of 
4.4  between  the  96  and  22  nm  samples  shortens  the  recovery  time  by  a 
factor  of  120  at  the  half-anneal  value  and  by  a  factor  of  260  at  the  three- 
quarter-anneal  point.  The  difference  in  these  results  is  caused  by  a  slight 
change  in  the  shape  of  the  recovery  curves  for  figure  5.  Specifically,  the 
recovery  curves  for  the  thicker  oxides  are  stretched  out  (more  dispersive) 
with  respect  to  the  curves  for  the  thinner  oxide.  Although  this  change  is 
relatively  small,  of  second  order  in  effect,  the  data  may  indicate  a  slight 
variation  in  the  a-value  with  thickness  for  the  as-grown  samples. 

As  noted  earlier,  at  low  temperature  and  field,  very  little  recovery  due  to 
hole  transport  occurs  on  the  time  scale  of  the  measurements.  However,  as 
may  be  discerned  from  the  lowest  temperature  recovery  curves  on  figures 
4,  5,  and  7  and  from  the  3-  and  4-MV/cm  curves  on  figure  9,  there  appears 
to  be  a  relatively  small  linear  (on  a  semilog  plot)  component  of  recovery 
that  constitutes  about  10  percent  of  the  recovery  at  8(H)  s,  and  which  seems 
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to  be  nearly  field  and  temperature  independent.  Hence,  a  different  mecha¬ 
nism  than  hole  transport  is  implied  for  this  component  of  recovery.  A 
possible  mechanism  is  that  of  removal  of  the  holes  generated  near  the  Si 
interface  by  electrons  tunneling  into  the  oxide  from  the  Si  and  annihilating 
these  holes.  This  is  made  more  plausible  by  the  79  K  recovery  data  1 17)  at 
1  MV/cm  shown  plotted  in  figure  12  for  the  three  oxide  thicknesses  of 
figure  10.  The  relative  fraction  of  recovery  due  to  the  log  t  component 
increases  as  the  oxide  thickness  decreases  in  just  about  the  proportion 
expected  for  a  constant  tunneling  distance  (at  a  given  time)  for  all  the 
samples.  Based  on  the  fraction  of  recovery  at  800  s,  we  estimate  that 
essentially  all  the  holes  within  ~3  nm  of  the  Si  interface  have  been  re¬ 
moved  by  tunneling  electrons  at  this  time.  Although  the  flatband  voltage 
recovery  due  to  the  tunneling/recombination  is  relatively  small -at  least  in 
the  thicker  samples-its  effect  still  has  to  be  accounted  for  in  some  of  the 
detailed  analysis  of  the  transport  component  to  be  performed  below. 


Figure  12.  Normalized 
low-icmpcrauire  (80  K) 
rial  hand  voltage  recovery 
data  for  etchcd-back  oxide 
samples  at  1 -MV/cm 
oxide  field,  indicating 
tunneling  recombination 
mechanism  for  recovery 
under  these  conditions. 
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5.  Analysis  and  Discussion  of  Results 


In  this  section  the  experimental  data  presented  in  the  previous  section  are 
analyzed  in  depth  to  establish  the  unusual  nature  of  the  delayed  hole 
transport  in  Si02.  After  first  examining  the  asymptotic  behavior  of  the 
charge  response  at  early  and  late  times,  we  quantitatively  characterize  the 
temperature,  field,  and  thickness  dependencies  of  the  response.  In  all 
cases  the  universality  of  the  time  response  with  respect  to  changes  in  the 
macroscopic  variables  (temperature,  field,  and  thickness)  is  presented  by 
replotting  the  data  in  terms  of  a  scaled  time  variable. 

5.1  Asymptotic  Analysis 

As  discussed  in  section  2  the  asymptotic  behavior  of  AVrB  provides  a 
critical  test  of  the  applicability  of  the  stochastic  transport  model.  If 
AVfb(0+)  is  the  initial  flatband  shift  immediately  after  the  radiation  pulse, 
then  the  relaxation  at  early  times  is  predicted  to  vary  as  IAV'r7?(0+)  - 
AVFB(t)\  ~  r  ,  whereas  at  long  time  A VFB(t)  is  predicted  to  decay  to  zero 
as  ra.  Hence,  if  the  CTRW  is  applicable,  a  log-log  plot  of  \AVFB(0+)  - 
AVFB(t)\  versus  t  at  early  time  should  yield  a  straight  line  of  slope  2a,  and 
a  log-log  plot  of  AVFB(t)  versus  t  at  late  time  should  yield  a  straight  line 
with  slope  -a. 

Figure  1 3  shows  the  results  of  applying  the  asymptotic  analysis  to  charge 
relaxation  data  taken  from  figures  4  to  7.  For  both  the  early-  and  late-time 
regimes,  three  sets  of  representative  data  were  analyzed  corresponding  to 
different  conditions  of  field  and  temperature,  which  are  noted  on  the 
figure.  In  all  cases  the  data  are  normalized  to  AVFB{Q+),  and  time  is  scaled 
in  units  of  ts  defined  as  follows:  In  the  early-time  regime,  ts  is  taken  as  the 
time  required  for  5-percent  recovery  (AVrB(,t)/AVFB(0+)  =  0.95),  and  in  the 
late-time  regime  the  time  is  scaled  to  the  time  for  95-percent  recovery 
(AVFB(t)/AVFB(0+)  =  0.05).  Hence,  the  universality  of  the  data  is  also 
examined  in  the  asymptotic  regimes. 

The  log-log  plot  of  the  early-time  response  is  shown  in  figure  13a.  The 
data  have  been  corrected  for  the  component  of  recovery  due  to  tunneling 
of  electrons  from  the  Si  substrate  to  trapped  holes  near  the  interface  (see 
fig.  12).  The  solid  line  is  a  least-squares  fit  to  the  data  between  5  x  10  2 
and  10  in  scaled  time  ts:  it  yields  a  slope  of  0.48  ±  0.04.  For  t/ts  >  10  the 
data  systematically  deviate  from  this  fit,  indicating  that  this  time  is 


Figure  13.  Asymptotic  analysis  of  flatband  voltage  recovery  data  for  96.5-nm  oxide  at 
several  values  of  field  (positive  bias)  and  temperatures,  (a)  Early-time  data  with  lime 
scaled  to  the  lime  at  which  5-percent  recovery  occurs,  (b)  Late-time  data  with  time 
scaled  to  the  95-percent  recovery  time.  Early-  and  lale-limc  slopes  are  consistent  with 
CTRW  model  with  dispersion  parameter  a  =  0.25  ±  0.03. 

approximately  the  upper  limit  to  the  validity  of  the  early-time  approxima¬ 
tion.  The  response  for  Eox  =  1  MV/cm  deviates  from  the  power  law  form 
approximately  a  half  decade  earlier  in  time  than  the  response  curves  for 
the  high  fields,  indicating  some  deviation  from  strict  universality  at  the 
lower  field.  This  may  be  attributed  either  to  effects  of  space-charge 
perturbation  on  the  oxide  field  or  to  some  intrinsic  change  in  the  nature  of 
the  hopping  paths  as  a  function  of  field. 

The  late-time  log-log  plots  are  shown  in  figure  13b.  The  experimental 
data  at  late  times  reveal  the  existence  of  a  remanent  charge,  corresponding 
to  about  1.8  percent  of  the  originally  generated  charge  still  present  hours 
after  irradiation,  which  evidently  represents  "permanent”  or  long-lived 
trapping  of  some  of  the  holes.  This  limiting  value  of  the  charge  relaxation 
was  subtracted  from  A VrB(t)  to  determine  the  late-time  behavior  of  those 
carriers  still  undergoing  stochastic  transport.  The  resulting  plot  of 
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logl(AV//7fl(r)  -  AV7fl)/AV/rfl(()+)l  versus  log  t  is  linear  over  at  least  four 
decades  in  time  with  a  slope  of  -0.25  ±  0.03.  Here  strict  universality 
(within  experimental  resolution)  is  obeyed  over  the  entire  time  range 
shown.  Departure  front  the  power  law  fomt  is  observed  for  tlts  <  3  x 
l(f2.  Now,  the  CTRW  stochastic  model  predicts  that  the  magnitude  of  the 
(negative)  slope  at  late  time  is  equal  to  the  (x  value  and  should  be  half  the 
value  of  the  slope  at  early  time.  Thus  the  relaxation  data  in  the  asymptotic- 
regimes  are  in  good  agreement  with  these  predictions  with  a  in  the  range 
from  0.22  to  0.27.  For  purposes  of  comparing  the  data  with  calculated 
response  curves  below,  we  take  a  =  0.25.  For  most  purposes,  variations  in 
a  of  ±0.03  or  so  are  relatively  insignificant. 

5.2  Temperature  Dependence 

We  first  examine  the  universality  of  the  hole  transport  with  respect  to 
temperature.  To  do  this  we  replot  in  figure  14  the  charge  relaxation  data 
of  figures  4  to  7  in  terms  of  a  scaled  time  variable,  which  we  take  to  be 
r1/2,  the  time  at  which  half  recovery  occurs.  The  solid  curve  in  each  case  is 
the  recovery  curve  calculated  from  approximate  solutions  of  the  CTRW 
model  for  a  =  0.25.  The  calculated  curve  is  also  scaled  to  the  half  recov¬ 
ery  time.  The  shape  of  the  calculated  curves  is  a  function  only  of  a;  the 
a-value  of  0.25  was  chosen  based  on  the  asymptotic  analysis  described  in 
the  previous  section.  For  the  case  of  E()X  =  -1  MV/cm,  the  solid  curve  is 
the  calculated  curve  for  negative  bias.  For  comparison,  the  predicted 
curve  for  positive  bias  is  also  shown  on  the  same  plot  as  the  dashed  curve, 
with  the  time  scaled  to  the  half  recovery  time  for  negative  bias.  Both  the 
shift  in  time  scale  of  -1-1/2  decades  and  the  difference  in  shape  of  the 
curves— primarily  evident  during  the  initial  recovery— result  from  the 
opposite  directions  of  charge  displacement  in  the  two  cases.  The  flatband 
voltage  reflects  the  first  moment  of  the  charge  density  relative  to  the  gate 
(metal)  electrode,  and  hence  the  recovery  curves  are  not  symmetric  with 
respect  to  bias  polarity. 

The  overall  universality  of  the  transport  data  with  temperature  is  evident. 
The  shapes  of  the  recovery  curves  are  essentially  the  same,  independent  of 
temperature;  only  the  time  scale  of  the  transport  changes  with  temperature. 
In  all  cases,  the  data  with  respect  to  the  scaled  time  trace  out  a  universal 
curve  over  10  or  more  decades  in  time,  with  data  points  in  several  in¬ 
stances  from  as  many  as  five  temperatures  overlapping  to  a  considerable 
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Figure  14  (cont’d).  Fixed  oxide  fields  of  (c)  l  MV/cm  and  (d)  5  MV/cm.  Solid  curves 
are  CTRW  model  calculations  for  a  =  0.25. 

degree  at  the  same  value  of  scaled  time.  There  appears  to  lie  a  tendency 
toward  increased  universality  with  increased  field,  as  discerned  from  a 
decrease  in  the  scatter  of  the  data  at  the  higher  fields.  This  may  be  due  to 
a  decrease  in  the  relative  perturbation  of  the  applied  field  due  to  space- 
charge  effects  as  the  field  is  increased.  Specifically,  based  on  the  actual 
values  of  the  Hatband  voltage  shifts,  for  an  applied  field  of  1  MV/cm  the 


actual  field  in  the  oxide  varies  initially  by  as  much  as  ±20  percent  because 
of  the  fields  generated  by  the  induced  hole  distribution.  At  the  higher 
fields  the  perturbation  is  much  smaller,  being  ~±5  percent  at  5  MV/cm.  It 
is  also  apparent  that  there  is  generally  very  good  agreement  between  the 
calculated  response  curve  and  the  data  over  the  scaled  time  regimes  shown 
(10  or  more  decades).  Only  at  Eox  =  1  MV/cm  is  there  a  marked  dis¬ 
crepancy  between  the  data  and  the  calculated  curve  for  early  scaled  time. 
This  is  most  likely  due  either  to  the  space-charge  perturbations  or  to  the 
tunneling  component  of  the  recovery  (fig.  12),  which  has  not  been  sub¬ 
tracted  from  the  data  here.  As  seen  from  figure  12,  the  tunneling  compo¬ 
nent  results  in  a  linear  increase  in  A VfR  with  log  time,  independent  of 

~  —4 

temperature,  which  is  clearly  evident  in  figure  14a  for  r/r]/2  <10  . 
Furthermore,  the  tunneling  component  ought  to  be  considerably  more 
pronounced  under  positive  bias  than  under  negative  bias,  because  under 
negative  bias  the  region  of  the  oxide  near  the  Si  interface  should  be 
essentially  devoid  of  holes  because  the  transport  of  holes  is  directed  away 
from  the  Si  interface.  This  is  borne  out  by  the  much  improved  agreement 
in  the  early-time  regime  between  the  calculated  response  curve  and  the 
data  for  Eox  =  -1  MV/cm  as  seen  in  figure  14b.  Here  there  is  no  discern¬ 
ible  tunneling  component. 

We  conclude  that,  except  for  the  tunneling  component  of  recovery  which 
could  be  accounted  for  quite  easily  by  use  of  figure  12,  the  theoretical 
curves  based  on  the  CTRW  with  a  single  value  of  a  provide  a  good 
description  of  the  data  over  the  entire  temperature  range  studied  and  for 
the  field  range  between  1  and  5  MV/cm.  Furthermore,  the  agreement 
exists  over  a  very  remarkable  range  of  scaled  time.  The  high  degree  of 
universality  over  the  entire  temperature  range  studied,  along  with  the 
highly  dispersive  nature  of  the  transport,  strongly  indicates  that  the 
stochastic  features  (the  dispersion)  of  the  transport  are  not  determined  by 
fluctuating  energy  levels  of  the  localized  hopping  sites,  but  rather  primar¬ 
ily  by  fluctuations  in  the  spatial  overlap  integrals  for  hopping.  We  con¬ 
clude  that  the  stochastic  parameter  a,  which  reflects  the  amount  of  disper¬ 
sion  in  the  transport,  is  essentially  temperature  independent  between  80  K 
and  room  temperature. 

It  may  be  instructive  to  compare  the  model  calculations  with  recovery  data 
in  real  time.  We  do  this  in  figure  15  for  an  oxide  field  of  3  MV/cm  (the 
data  of  fig,  6),  by  simply  shifting  the  time  scale  for  the  model  curve  to  the 
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Figure  1 5.  Comparison  of 
CTRW  model  calculations 
of  Hatband  voltage 
recovery  lor  dispersion 
parameter  value  of  0.25 
with  measured  response  in 
real  time  for  oxide  field  of 
3  MV/cm. 


real  time  appropriate  for  each  temperature.  Again,  the  agreement  between 
the  data  and  model  curves  is  striking. 

We  now  analyze  the  temperature  dependence  of  the  transport  quantita¬ 
tively  by  examining  the  variation  in  hole  transit  time  with  temperature  and 
determining  activation  energies  for  the  transport  process.  Figure  16 
presents  the  activation  energy  analysis  for  the  transport  data  of  figures  4  to 
7.  These  plots  are  generated  by  the  method  of  "cross-cuts."  We  assume 
that  for  each  applied  field  there  is  a  one-to-one  correspondence  between 
the  normalized  flatband  voltage  and  a  single  intrinsic  scaled  time,  inde¬ 
pendent  of  temperature,  as  the  universal  plots  would  indicate.  The  intrin¬ 
sic  scaled  time  is  simply  an  indicator  of  a  certain  amount  of  transport 
having  taken  place.  Then  we  plot  the  real  time  necessary  to  obtain  a  fixed 
value  of  A VrD  versus  1  fT.  Figure  16  shows  such  Arrhenius  plots  for  a 
series  of  normalized  A!  rB  values  from  a  value  at  which  little  transport  has 
occurred  (AV rB/AV  rp(0* )  =  0.90,  corresponding  to  early  scaled  time),  to  a 
value  where  most  of  the  transport  has  occurred  (A\  r/?/AV  F/?(0+)  =  0.10, 
corresponding  to  late  scaled  time).  Each  set  of  curves  is  for  a  fixed  value 
of  the  applied  field.  The  data  in  figure  16  were  corrected  for  the  tunneling 
component  of  the  recovery,  a  correction  of  practical  significance  only  at 
low  temperature. 

The  first  important  point  to  be  noted  from  figure  16  is  that  there  are  two 
distinct  temperature  regimes:  a  high-temperature  regime  in  which  the 
transport  is  distinctly  activated  and  a  low-temperature  regime  in  which  the 
transport  becomes  essentially  not  thermally  activated  (or  activated  with  an 


energy  only  of  the  order  of  0.01  eV).  The  transition  from  activated 
behavior  to  nonactivated  behavior  seems  to  begin  at  about  160  K  going 
down  in  temperature,  and  based  upon  the  high  field  data  (fig.  16d),  seems 
to  be  complete  around  1  (X)  K.  The  transition  is  barely  observed  to  be 
beginning  at  low  field  (±1  MV/cm)  because  even  for  only  10-percent 
recovery  to  occur  at  140  K  requires  ~  1 0  s,  which  is  the  upper  limit  of  our 
measurement  time.  As  the  field  increases  the  transport  speeds  up  rapidly, 
and  for  Eox  =  5  MV/cni  both  temperature  regimes  are  plainly  encompassed 
by  the  measurements  even  for  75-percent  recovery  (AVrD/AVrD(0+)  = 
0.25).  The  3-MV/em  data,  of  course,  represent  an  intemiediate  case. 

The  temperature  dependence  of  the  transport  observed  here  is  a  classic 
signature  of  polaron  hopping  [40-451,  as  discussed  in  an  earlier  paper  1 1 1  j 
and  in  section  2.5.  Most  polaron  theories  predict  a  changeover  from 
thermally  activated  transport  to  nonactivated  transport,  or  to  a  process 
characterized  by  a  much  lower  activation  energy  corresponding  to  the 
emission  or  absorption  of  a  single  phonon,  as  the  temperature  is  lowered 
below  about  1/3  or  1/4  of  the  Debye  temperature.  This  is  consistent  with 
our  results  as  ~  6(X)  K  for  Si02.  Other  mechanisms  such  as  variable 
range  hopping  cannot  be  completely  ruled  out  based  on  the  temperature 
data  alone.  However,  further  arguments  for  a  small  polaron  transport 
process  come  from  recalling  the  high  degree  of  universality  of  the 
transport  with  temperature:  also,  the  CTRW  model  curve  with  a  single 
value  of  the  stochastic  a-parameter  provides  a  good  description  of  the  data 
over  more  than  10  decades  in  time,  including  both  the  early-time  (low- 
temperature)  and  late-time  (high-temperature)  regimes.  That  such  univer¬ 
sality  is  exhibited,  in  spite  of  the  fact  that  the  activation  energies  for  the 
transport  are  very  different  at  low  and  high  temperatures,  strongly  suggests 
that  the  hopping  sites  involved  in  the  hole  transport  are  the  same  at  all 
temperatures  and  times.  Only  the  microscopic  details  of  the  temperature- 
activation  process  by  which  a  hole  makes  the  transition  between  two  sites 
vary  with  temperature.  The  stochastic  features  of  the  transport,  on  the 
other  hand,  are  determined  by  the  fluctuations  in  the  spatial  overlap 
integrals  for  hopping.  This  is  all  consistent  with  hole  transport  via  small 
polaron  hopping  between  localized  sites. 

Focussing  attention  now  on  the  high-temperature  activated  regime,  several 
points  should  be  noted.  First,  during  the  transport  time  for  10  to  75 
percent  of  the  recovery  to  occur,  the  activation  energies  for  each  field 


value  as  ascertained  from  the  slopes  of  the  Arrhenius  plots  have  a  fairly 
narrow  spread  in  values,  of  the  order  of  0.03  eV.  Only  at  late  scaled 
time— at  90-percent  recoverv--is  there  a  significant  change  in  the  activation 
energy,  an  increase  in  all  cases  of  -0.1  eV.  This  would  indicate  that  for 
most  of  the  recovery  process  the  hole  transport  reflects  a  single  activation 
energy,  or  more  generally,  the  holes  are  in  equilibrium  over  a  distribution 
of  activation  energies.  This  would  tend  to  rule  against  a  conventional 
multiple-trapping  model  for  most  of  the  transport,  as  in  this  case  the 
activation  energy  should  increase  progressively  with  the  fraction  of 
transport,  or  scaled  time.  The  increase  in  the  activation  energy  observed  at 
late  time  would  indicate  that  deeper  traps  are  coming  into  play  at  this 
stage.  Hence,  for  the  later  stages  of  the  transport  a  distribution  in  activa¬ 
tion  energies  may  be  beginning  to  significantly  affect  the  transport. 
Another  possible  explanation  for  the  change  in  A  at  late  time  is  that  the 
hole  distribution  at  late  time  is  skewed  toward  the  negative  interface  and, 
hence,  may  be  reflecting  a  change  in  the  trap  energy  distribution  near  an 
interface.  Note,  however,  that  under  both  positive  and  negative  polarities 
(±1  MV/cm)  the  values  and  changes  in  A  with  time  are  essentially  the 
same;  therefore,  the  energy  distributions  would  have  to  be  very  similar 
near  both  interfaces.  This  would  appear  somewhat  fortuitous,  and  sug¬ 
gests  to  us  that  the  more  likely  explanation  is  that  the  same  distribution  of 
activation  energies  is  involved  throughout  the  transport.  Initially,  when 
the  number  of  holes  in  the  sample  is  relatively  large,  they  rapidly  come 
into  equilibrium  with  the  distribution  of  activation  energies.  The  number 
of  holes  interacting  with  deeper  energy  levels  remains  nearly  constant  (i.e.. 
these  levels  are  saturated)  and  constitutes  a  small  fraction  of  the  total  hole 
population.  But  as  the  transport  progresses  until  most  of  the  holes  have 
been  collected  at  the  cathode,  the  number  of  hops  reflecting  the  larger 
activation  energies  becomes  an  increasingly  larger  fraction  of  the  total 
number  of  holes,  which  shows  up  as  an  increase  in  the  observed  average 
activation  energy.  This  is  certainly  a  manifestation  of  multiple  trapping, 
and  the  above  arguments  may  offer  a  basis  for  reconciling  seemingly 
conflicting  observations  reported  previously 

The  next  point  to  be  noted  about  the  high  temperature  activated  regime  of 
figure  16  is  the  field  dependence  of  the  activation  energy  There  is  a 
marked  decrease  in  A  with  increasing  electric  field  This  is  shown  more 
clearly  in  figure  17  where  we  plot  the  activation  energies  obtained  from 
figure  16  at  various  scaled  limes  versus  field.  The  marked  change  in  the 


value  of  A  during  the  latter  stage  of  the  transport  process  is  clearly  shown. 
For  each  field  value  the  points  for  10-,  25-.  50-,  and  75-percent  recovery 
all  cluster  together  within  experimental  accuracy,  whereas  the  A-values  at 
90-percent  recovery  (A\'rB/AVrB(Q+)  =  0.10)  all  lie  about  0.1  eV  higher. 
In  both  cases,  however,  there  is  a  linear  decrease  in  A  with  increasing  field 
with  a  slope  of  about  -0.05  eV/MV/cm.  More  precisely,  based  on  the 
straight-line  fits  to  the  points  in  figure  17,  the  field  dependence  of  A  can  be 
expressed  by 

A(Em)  =  A-hE,„  .  (7) 


Figure  17.  Activation 
energies  obtained  from 
figure  16  for  several 
fractions  of  recovery 
plotted  versus  oxide  field. 


where  h  -  0.05  eV/MV/cm  and  Eox  is  in  megavolts  per  centimeter.  The 
zero  field  activation  energy  Atl  =  0.64  eV  over  most  of  the  transport  time, 
up  through  AVrB/A\'rB[()+)  =  0.25,  but  changes  to  At/  =  0.75  eV  during  the 
final  transport  stages.  The  linear  decrease  in  A  with  field  is  consistent  with 
a  hopping  transport  process  where  b  =  ea/2,  with  a  =  1  nm  being  the 
average  hopping  distance  in  the  forward  (field)  direction.  We  now  con¬ 
sider  the  field  dependence  of  the  transport  more  explicitly. 
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Field  Dependence 

The  field  dependence  of  the  hole  transit  time  and  its  interpretation  in  terms 
of  polaron  hopping  was  discussed  in  an  earlier  paper  1 15|.  However,  there 
the  universality  of  the  transport  data  was  not  shown  explicitly.  Here  we 
first  focus  on  the  universality  feature  and  then  briefly  review  the  analysis 
of  the  dependence  of  hole  transit  time  with  field. 

The  universality  of  the  transport  with  field  at  T  =  79  K  is  shown  in  figure 

18.  There,  in  the  same  manner  as  done  for  the  temperature-dependent 

data,  the  data  of  figure  9  are  replotted  in  terms  of  a  time  scaled  to  r)/2,  the 

time  required  for  half  recovery  of  AV^g.  Again  the  recovery  curves  very 

nearly  trace  out  a  single  "universal"  curve  in  terms  of  the  scaled  time 

variable  over  essentially  the  entire  time  of  transport.  Some  amount  of 

sharpening  of  the  recovery  curve  is  apparent  at  6  MV/cm  relative  to  the 

lower  field  data  (in  the  range  tltx/7  =10  to  10°),  indicating  that  perhaps 

another  effect  or  process  may  be  coming  into  play  at  very  high  field,  but 

this  seems  to  be  a  minor  effect  overall.  The  solid  curve  is  the  same  model 

curve  for  a  =  0.25  (and  for  positive  bias)  as  used  in  figure  14.  Except  for 

—  -2 

a  relatively  small  discrepancy  at  early  scaled  time  (<10  ),  which  again 
may  be  attributed  to  the  tunneling  component  of  recovery,  the  agreement 
between  the  data  and  the  calculated  curve  is  very  good. 


CTRW  MODEL 

(*  =  0  25)  ^ _ Pc 


O  5  MV/cm 
□  4  MV/  cm 

0  3  MV/  cm 


SCALED  TIME  (I  1 ,  2) 

Figure  18.  Normalized  llalhand  voltage  recovery  data  of  figure  9  replollcd  versus  scaled  lime 
indicating  universality  ol  recovery  with  field  at  constant  temperature  (79  K)  and 
thickness  (96.5  mil ».  Solid  curve  is  (TRW  nuxlcl  calculation  lor  a  =  0.25. 


In  earlier  work  [  I5|.  the  field-dependent  data  (fig.  X  and  9)  were  analyzed 
in  terms  of  the  semiclassical  small  polaron  hopping  mobility: 


|i(£(<)  =  p<0)  sinh([k’E  >ta/2)/(PeE  i(<i/2>  (Xi 

where  Eox  is  the  oxide  electric  field,  p  =  (A7>  '.  a  is  an  average  hopping 
distance,  and  |i(0)  is  the  low  field  mobility.  We  take  the  half-recovery 
time  t]f2  as  a  measure  of  the  transit  time  and  hence  as  a  measure  of  the 
average  inverse  mobility.  (Because  of  the  anomalous  nature  of  the  disper¬ 
sive  hole  transport  in  SiO,.  the  use  of  the  term  'mobility  has  limited 
usefulness  in  the  conventional  sense,  since,  as  discussed  in  sect.  2.  it  turns 
out  that  the  effective  average  mobility  is  both  time  and  thickness  depend 
ent.  We  are  using  the  concept  in  a  more  general  sense  here:  more  specifi 
tally,  we  are  assuming  the  hole  transit  time  simply  reflects  the  field 
dependence  of  the  inverse  average  intersite  hopping  rate.)  For  the  case  ot 
large  argument  (|leE  <f/2  »  I ).  which  is  true  for  the  fields  and  tempera¬ 
tures  here,  equation  (X)  reduces  to  the  exponential  form 

/,/,<£ )t>  =  f,n(0)exp(-PeE„£//2)  .  (Ml 

where  r,/2(0)  is  a  constant.  In  figure  19  we  show  log  t]  . ,  plotted  versus  F.  t 
with  temperature  as  a  parameter.  The  data  have  been  corrected  for  the 
tunneling  component.  The  exponential  dependence  of  the  recovery  with 
field  is  evident,  with  the  magnitude  of  the  slope  increasing  as  the  tempera¬ 
ture  is  lowered.  If  the  semiclassical  hopping  formula  is  applicable,  the 
field  slope  should  be  ca/lkT. 

Figure  20  shows  the  field  slopes  plotted  versus  1/7.  The  circles  are  the 
values  of  the  slopes  taken  from  figure  19  for  the  times  required  for  half 
recovery.  The  straight  line  is  a  fit  through  these  points  and  the  origin. 
According  to  the  hopping  mobility  formula,  the  slope  of  the  line  yields  a 
value  of  0.9  nm  for  the  average  hopping  distance.  Also  shown  in  figure  20 
are  points  taken  from  other  plots  similar  to  figure  19  where  other  fixed 
fractions  of  charge  recovery  (12.5.  25.  and  75  percent)  were  used  as  a 
measure  of  the  hole  transport  speed.  Only  transit  time  versus  field  plots 
with  at  least  three  field  points  are  included.  Though  there  is  some  scatter 
in  these  points,  all  the  data  are  consistent  within  experimental  resolution 
with  the  straight  line  through  the  half  recovery  points.  From  the  uncer- 
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tainties  in  the  data  and  the  analysis,  we  estimate  error  limits  on  the  average 
hopping  distance  of  ±0.2  nm.  We  note  that  an  average  hopping  distance  of 
0.9  nm  is  consistent  with  the  reduction  in  activation  energy  we  have 
observed  with  increasing  field  in  the  high-temperature  activated  hopping 
regime  (see  fig.  17). 


5.4  Thickness  Dependence 

For  completeness  we  briefly  review  the  analysis  of  the  thickness  depend¬ 
ence  of  the  hole  transport  discussed  elsewhere  [17|.  First  recall  that 
because  of  the  anomalous  nature  of  the  hole  transport  process,  the  CTRW 
model  predicts  that  the  dependence  of  the  transit  time  on  oxide  thickness 
Lox  is  superlinear,  of  the  form  L  'j“.  In  the  wet  oxides  used  in  the  present 
investigation  we  have  found  that  a  =  0.25  ±  0.03  from  the  asymptotic 
analysis.  Further,  the  calculated  model  curves  based  on  this  value  of  a  are 
in  excellent  agreement  with  the  recovery  data  over  a  wide  range  of 
temperature  and  field  values  and  over  many  decades  in  time.  Therefore, 
based  on  the  model,  the  hole  transit  time  through  the  oxide  should  vary 
approximately  as  L*x . 

To  compare  the  theory  and  experiments  on  a  quantitative  basis,  we  show 
in  figure  21  log-log  plots  of  the  recovery  time  versus  oxide  thickness  for 


the  recovery  data  of  figures  10  and  11,  for  the  etched-back  and  as-grown 
sets  of  samples,  respectively.  The  relative  displacement  between  the  two 
plots  is  not  significant,  since  the  data  for  the  two  plots  were  obtained  under 
different  conditions.  For  the  etched-back  samples,  the  data  were  recorded 
at  79  K  and  5  MV/cm  (fig.  10)  and  the  time  for  A\'rB  to  reach  half 
recovery  was  used  as  a  measure  of  the  transit  time.  For  the  as-grown 
samples,  the  data  were  taken  at  220  K  and  1  MV/cm  (fig.  1 1),  and  the  time 
for  AVrB  to  reach  75-percent  recovery  was  used  as  a  measure  of  the  transit 
time.  (We  used  the  75-percent  recovery  point  in  this  latter  case  in  order  to 
have  a  reliable  horizontal  cut  through  all  the  recovery  curves  in  fig.  11.) 
The  importance  of  figure  21  is  that  it  shows  that  the  log  recovery  time 
versus  log  thickness  data  can  be  fit  reasonably  well  by  straight  lines  in 
both  cases,  with  slopes  of  4.2  and  4.0  for  the  etched-back  and  as-grown 
cases,  respectively.  These  values  agree  very  well  with  the  prediction  of  an 
L*  dependence  for  this  oxide. 

Next  we  show  the  universality  of  the  data  with  respect  to  thickness.  In  the 
usual  manner  we  replot  in  figure  22  the  recovery  data  for  figure  10  for  the 
etched-back  samples  in  terms  of  the  scaled  time  at  which  half  recovery 
occurs.  Again  the  solid  curve  is  the  calculated  recovery  curve  for  a  = 
0.25.  The  universal  feature  of  the  data  is  evident,  and  the  calculated  curve 
is  in  good  agreement  with  the  data  over  ~8  decades  in  time.  Essentially 
the  same  result  holds  true  for  the  as-grown  data. 
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Figure  2 1 .  Log  of  hole  transit  lime  plotted  versus  log 
of  oxide  thickness.  For  the  clched-back  samples,  the 
time  to  achieve  half-recovery  was  used  as  a  measure 
of  transit  lime.  For  the  as-grown  samples,  75-percent 
recovery  time  was  used  because  of  scarcity  of  data 
near  half-recovery  for  the  2 1 ,8-nm  oxide. 


Scaled  lime  (t/l 

Figure  22.  Normalized  Hatband  voltage  recovery  data  ol 
figure  10  (clched-back  oxide  samples)  plotted  against 
scaled  lime  (i/t]f2).  indicating  universality  with  oxide 
thickness  at  constant  temperature  (80  K)  and  field 
(5  MV/cm).  Solid  curve  is  calculated  recovery  curve  for 
a  =  0.25. 


6.  Summary 


In  this  report  we  have  demonstrated  that  the  hole  transport  in  SiO,  exhibits 
universality  in  terms  of  a  scaled  time  variable  over  wide  ranges  of 
temperature,  field,  and  thickness  lurthermore.  the  C'l'RW  model  with  a 
single  value  of  0.25  for  the  stochastic  tx  parameter  provides  a  good 
description  of  the  recovery  data  over  a  remarkable  range  (-10  decades)  in 
scaled  time  for  the  range  of  temperature,  field,  and  thickness  investigated 
here.  This  is  true  at  least  for  the  clean,  hardened  oxides  used  in  these  in¬ 
vestigations,  i.e.,  oxides  which  have  little  (<2  percent)  long-term  deep  hole 
trapping.  This  is  a  crucial  result  which  indicates  rather  strict  universality 
of  the  hole  transport  with  respect  to  temperature,  field,  and  thickness  in 
clean  SiO,.  Changing  these  parameters  simply  changes  the  time  scale  for 
the  response  but  not  the  shape  or  overall  dispersion  of  the  response  curves. 
Summarizing  in  more  quantitative  terms,  a  response  function  (radiation- 
induced  current  or  voltage  shift)  may  be  characterized  by  the  functional 
form  F(a:t/tJ,  where  the  disorder  parameter  a  describes  the  dispersion  tor 
shape)  of  F  and  is  independent  of  /  ,  F  .  and  L  .  These  parameters  enter 
only  in  the  characteristic  time  scale.  /?.  for  the  response:  fv  =  /t(7\  Etn. 
Lox).  The  universality,  particularly  with  respect  to  temperature,  has  im¬ 
plications  for  the  nature  of  the  transport,  namely,  that  the  transport 
proceeds  via  hopping  between  localized  sites.  The  dispersion  originates 
primarily  from  a  distribution  of  intersite  transfer  integrals,  as  would  arise 
either  from  a  random  spatial  distribution  of  hopping  sites,  or  perhaps  from 
a  distribution  of  bond  angles  in  the  SiO,  network.  However,  the  field  de¬ 
pendence  of  the  transport-  specifically  the  field  dependence  of  the  activa¬ 
tion  energy  -points  to  a  hopping  prx>cess  with  an  average  hopping  distance 
in  the  field  direction  of  ~1  nm.  The  actual  intersite  transfer  process  seems 
most  likely  to  be  small  polaron  hopping  as  indicated  by  the  transition  from 
thermally  activated  transport  above  -140  K  to  an  essentially  nonactivated 
transport  at  lower  temperature. 

The  characteristic  time  scale  for  the  transport,  which  is  directly  related  to 
the  distribution  of  transit  times  of  the  holes  across  the  oxide  film,  has  been 
found  to  be  expressed  in  terms  of  temperature,  field,  and  thickness  as 


where  a  is  the  average  hopping  distance,  t"s  is  a  constant,  the  field 
dependent  activation  energy  is 


me, J  =  a 


W  Eox 


(ID 


and  A()  is  the  low'  field  limit  of  A (Eox).  For  the  particular  oxides  studied 
here,  Ao  =  0.64  eV  and  a  =  0.9  nm  which  leads  to  about  a  0.05-eV  reduc¬ 
tion  in  A  per  MV/cm  increase  in  field.  Also  for  this  oxide,  a  =  0.25  and, 
using  the  half-recovery  point  for  AVFB  as  a  measure  of  / ,  we  find  t”  =  3x 
l(r23  s.  With  these  parameter  values,  the  time  for  half-recovery  to  occur 
according  to  equation  ( 10)  for  a  100-nm  oxide  at  room  temperature  and  for 
a  1 -MV/cm  field  is  ~1  x  10  s  (see  fig.  4).  Finally,  we  note  that  in  the 
low-temperature  nonactivated  regime.  kT  in  equation  (10)  is  replaced  by 
vkTD  where  is  the  Debye  temperature  and  v  is  a  numerical  factor  on  the 
order  of  1/4.  Even  though  the  transport  is  not  thermally  activated  in  this 
temperature  regime,  it  is  still  field  activated,  with  rv  oc  exp  {-nqEj2vkTn). 
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